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The observation of the top (t) quark served as one of the major confirma- 
tions of the validity of the standard model (SM) of particle interactions |1I2| . 
Through radiative corrections of the SM, the mass of the top quark, along 
with that of the W boson [jy* , provide the best indication for the value of the 
mass of the hypothesized Higgs boson [3]. The mass of the W is known to a 
precision of < 0.1%, while the uncertainty on the mass of the top quark is at 
the 4% level Improvements in both measurements are required to limit 
the range of mass that the Higgs boson can assume in the SM, and, of course, 
to check whether that agrees with expectation. It is therefore important to 
develop techniques for extracting the mass of the top quark that can provide 
the sharpest values possible. 

Measurement of M t 

We report on a new preliminary measurement of the mass of the top quark 
from ti data in lepton+jets channels accumulated by the D0 experiment in 
Run-1 of the Tevatron. The luminosity corresponds to 125 events/pb, and 
this analysis is based on the same data sample that was used to extract the 
mass of the top quark in our previous publication 5 . Information pertaining 
to the detector and to the older analysis can be found in Refs. |Bj and [S], 
respectively. 

After offline selections on lepton transverse energies (Et > 20 GeV) and 
angles flr^ < 1.7, and \r] e \ < 2.0, for muon and electron channels, respec- 
tively), on jet transverse energies (Et > 15 GeV) and angles (I77I < 2.0), 
imbalance in transverse momentum (missing- Et > 20 GeV), and after ap- 
plying several less important criteria [jy, the event sample consisted of 91 
events with one isolated lepton and four or more jets. (Unlike the previous 
analysis, we do not distinguish between events that have or lack a muon as- 
sociated with one of the jets, signifying the possible presence of a &-quark jet 
in the final state.) The new analysis involves a comparison of the data with 
a leading-order matrix element for the production and decay process, and to 
minimize the effect of higher-order corrections we therefore restrict the study 
to events containing only four jets. This requirement reduces data sample to 
71 events. 

In the previous analysis, the four jets with highest Et were assumed 
to represent the four quarks in the events (two b quarks from the initial 
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Fig. 1. Discriminator D = P sig „ a i/(P S i g nai + Pbackg.) for the new analysis, and 
(inset) for the two versions of the analysis included in the previous D0 publication 




t ^ W + b decays, and the q'q decays of one of the W bosons). These, along 
with the lepton and the missing v (from the decay of the other W boson) 
were fitted to the kinematic hypothesis pp — > tt — > WWbb, subject to the 
constraints of overall momentum-energy conservation, the known mass of the 
W boson, and the fact that the unknown mass of the top quark (M t ) had 
to be identical for the top and antitop quarks in the event. With 12 ways 
to permute the jets, there were 12 possible fits, and the solution with the 
lowest x 2 (but < 10) was chosen as the best hypothesis, thereby defining 
the fitted mass mfu for the event, as well as the longitudinal momentum of 
the v. The same procedure was used for generating templates in variables 
of interest, as a function of input mass of the top quark. This was based on 
the Monte Carlo (MC) HERWIG program 0, which was used to generate 
events that were passed through full event-reconstruction in the detector 
The background, which consisted mainly of all-jet production (20%) and 
W+jets (80%) was also processed in a similar manner. The background from 
all-jet production was based on data, and the background from W+jets was 
based on events generated with VECBOS 0. A discriminant was formed 
based on differences in distributions in four variables predicted for tt signal 
and background, and this defined the probability that any event represented 
signal as opposed to background. A probability density was defined as a 
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function of the discriminant (D) and mja, and a comparison of data and 
MC via a likelihood was used to determine the most likely mass of the top 
quark. 

The new analysis is also based on a likelihood, but this likelihood is a 
function of all measured variables in the event, with the exception of the un- 
clustered energy in the calorimeter that is used to define the missing Et- Our 
method is similar to that suggested for ti dilepton decay channels JU]j an d 
used in previous mass analyses of dilepton events We define a probability 
for background purely in terms of the matrix element contained in VECBOS, 
and for signal in terms of the leading-order matrix element for tt production 
and decay. This is convoluted with a transfer function that relates objects at 
the parton level to the objects observed (fully reconstructed) in the detector. 
We assume that the observed electrons correspond to the produced electrons, 
and that the muons are smeared with their known resolution. The angles of 
the jets are assumed to reflect the angles of the partons in the final state, 
and we ignore any transverse momentum for the incident partons. We use a 
large MC sample of ti events (generated with masses between 140—200 GeV 
in HERWIG, and processed through the D0 detector-simulation package) to 
determine a set of ten parameters that correlate any jet Et with its parton 
value. The parameters used for b quarks are different than for the lighter 
quarks, and there are therefore 20 parameters in all. 

We can write the probability for ti production in terms of the following 
elements: 



where the sum is over all 12 permutations of the jets, and all possible values 
of v momenta. |M| 2 is the matrix element for the process, f(qi) and J(q2) 
are the parton distribution function (PDF) for the incident particles, is 
the phase-space factor for the 6-object final state and W(x, y) correspond to 
a function that parameterizes the mapping between parton level quantities 
y and detector measurements x. With two incident parton energies (we take 
these partons to be quarks, and ignore the ~ 10% contribution from gluon 
fusion), and six objects in the final state, the integrations over the essentially 
fifteen sharp variables (3 components of lepton momentum, eight jet angles, 
and four <5- functions representing energy- momentum conservation), leaves five 
integrals that must be performed to obtain the probability that any event 
represents ti production for some specified value of top mass m t . Four of 
the variables chosen for the remaining integrations, namely the masses of the 
W bosons and of the top quarks in the event, are economical in CPU time, 
because the value of \M\ 2 is essentially negligible except at the peaks of the 
four Breit-Wigner terms in the matrix element. 
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Fig. 2. Distribution for background probability calculated for the 71 4-jets ti can- 
didates (data points). Only the events to the left of the vertical line (\n[Pb ac kg] < 
— 11,0) are considered for further analysis. The data is compared with the expected 
results from MC-simulated samples (dashed histrogram) formed with a mixture of 
16 signal (left-hatched histrogram) and 55 background (right-hatched histogram) 
events. 

Introducing an analogous expression for the background, the likelihood 
as a function of mt can be written as: 



where A(x) is the acceptance of the detector in terms of the parton variables 
(and contains the transfer function), and the sum over the probabilities for 
the individual events reflects the product of the probabilities for final sample 
of N of events. The best values of a, representing m t , and the parameters Cj, 
are all defined by the most probable value of the likelihood. 

The main difference between this method and the previous analysis is 
that each event now has its individual probability as a function of the mass 
parameter. This probability, reflecting both signal and background, depends 
on all measured variables in the event (excepting unclustered energy), with 
well-measured events contributing more sharply to the extraction of the mass 
of the top quark than those poorly measured. The fact that so much of the 






Measuring M t in tt events. 



5 



information in each event is retained in the probability provides a far better 
separation of signal and background. Figure^displays the effective difference 
in the discriminating power for the parameter D in the current and in the 
previous analyses. In all these analyses, this parameter is not used to select 
regions of high signal-probability, but the comparison is instructive in showing 
the greater power of the present method. 




Fig. 3. a) Negative of the log of the likelihood as a function of the mass of the 
top quark for the 22 tt candidates in our final sample, b) Probability distribution 
determined from the likelihood, with the hatched area corresponding to the 68% 
probability interval, from which we determine M t = 179.9 ± 3.6 GeV. A 0.5 GeV 
upward bias correction has been included based on the results of MC ensemble 
tests. 



Studies on samples of HERWIG MC events used in the previous analy- 
sis indicate that the expected uncertainty using the new method would be 
capable of yielding almost a factor of two reduction in the uncertainty of 
the extracted top mass. However, the method appears to be sensitive to a 
systematic shift of the mass (of about 2 GeV) as a function of the a pri- 
ori sample purity, whose true value cannot be determined in the analysis. 
To minimize the impact of this limiting systematic uncertainty, we proceed 
as follows. Figure |2 shows the probability that each event corresponds to 
background (VECBOS), and we compare that to a sample of 71 MC events 
corresponding to the mixture of 16 signal and 55 background events found 
in the data sample. From MC studies we fibd that eliminating from further 
consideration in the analysis those events that have a high probability of 
being background, reduces markedly the dependence on the extracted mass 
on sample purity. We consequently retain for further study only those events 
that have a poor probability of being background. This selection is shown in 
Fig and its imposition leaves a sample of only 22 events, 12 of which are 
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signal and 10 background, after the recalculation of the new minimum in the 
likelihood, the result of which is shown in Fig. [3] 

Another aspect of this analysis is that it provides a natural method for 
measuring any other parameter in the ti differential cross section. As an ex- 
ample we, show in Fig. 0] the likelihood as a function of the mass of the W 
boson (M\y) for the same sample of events. This suggests the very interest- 
ing possibility of checking the jet energy scale (JES) in this events, which 
corresponds to the largest contribution to the systematic uncertainty in the 
current and previous analyses. The next step in this analysis is, in fact, to 
utilize the information on JES contained in these events. However, in this pre- 
liminary result, we only estimate this uncertainty using the previous analysis 
0, where the JES calibration from 7 + jet events was used to match the en- 
ergy scale in the experiment to that in the MC simulation. Here we estimate 
the uncertainty in M t coming from JES as follows. The analysis is performed 
before and after applying this correction, and the difference defined as the 
systematic uncertainty due to JES. Table H summarizes all the systematic 
uncertainties in this preliminary result. 
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Fig. 4. a) Negative of the log of the likelihood as a function of the mass of the 
W boson for the 22 ti candidates in our final sample, b) Probability distribution 
determined from the likelihood, with the hatched area corresponding to the 68%. 
probability interval. 



Conclusion 

We have presented a new preliminary measurement of the mass of the top 
quark using a method that compares each individual event with the differ- 
ential cross section for ti production and decay. We obtain a significant im- 
provement over the statistical uncertainty of the previous measurement [5], 
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that is equivalent to having a factor of 2.4 more data. The new preliminary 
result is: 

M t = 179.9 ± 3.6(stat) ± 6.0(sys) GeV (3) 

The possibility of checking the value of the W mass in the same events 
offers the possibility of redoing the largest systematic uncertainty, of the same 
events provides a new handle on controlling the largest systematic, namely 
the jet energy scale. 

Table 1. Table of systematic uncertainties for the measurement of M t 



Model for tt 


1.5 


GeV 


Model for backgound(tT+jets) 


1.0 


GeV 


Noise and multiple interactions 


1.3 


GeV 


Jet energy scale 


5.6 


GeV 


Parton distribution function 


0.2 


GeV 


Acceptance correction 


0.5 


GeV 



Total 6.0 GeV 
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